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We present a study of the nucleation mechanism that allows the decay of the 

•rH . 

| metastable phase (trans-cisoid) to the stable phase (cis-transoid) in quasi 



one-dimensional non-degenerate polymers within the continuum electron- 
phonon model. The electron-phonon configurations that lead to the decay, 
i.e. the critical droplets (or transition state), are identified as polarons of 
the metastable phase. We obtain an estimate for the decay rate via thermal 
activation within a range of parameters consistent with experimental values 
for the gap of the cis-configuration. It is pointed out that, upon doping, the 
activation barriers of the excited states are quite smaller and the decay rate 
is greatly enhanced. Typical activation energies for electron or hole polarons 
are ~ 0.1 eV and the typical size for a critical droplet (polaron) is about 
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20A. Decay via quantum nucleation is also studied and it is found that the 

crossover temperature between quantum nucleation and thermal activation is 

of order T c < A0°K. Metastable configurations of non-degenerate polymers 

may provide examples for mesoscopic quantum tunneling, 
pacs No:71.10.+x;71.20.Hk;64.60.My 
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I. INTRODUCTION AND MOTIVATION 



Quasi-one dimensional conductors provide a fascinating wealth of physical phenomena 
that stimulated considerable theoretical and experimental study of these systems. The 
proposal [ffl that solitons play a central role in the electronic and transport properties of the 
degenerate polymer trans-(CH) x polyacetylene was later confirmed with measurements of 
the optical absorption in the doped system . Non-topological excitations corresponding 
to local lattice deformations that bind electron states, i.e. polarons, had also been recognized 
in these degenerate polymers 

Although the presence of topological solitons is associated with the ground state de- 
generacy of irans-polyacetylene, polarons are not a consequence of this degeneracy and are 
predicted to occur upon doping also in non-degenerate isomers ||. 

Non-degenerate isomers like cis-(CH) x were studied || by using a continuum electron- 
phonon model ||. These non-degenerate isomers feature a globally stable (lower energy) 
configuration (cis-transoid) and a locally stable but globally metastable (higher energy) 
configuration (trans-cisoid) . 

The properties of the polaron excitations of the stable (cis) phase were thoroughly studied 
within the continuum electron-phonon model and show the property of "confinement" , that 



is their energy grows linearly with the spatial size of the lattice distortion |p||T|,P|,|TU|] . Because 
of the lack of degeneracy, there are no stable topological soliton excitations in this case. 

The motivation for our study is to understand the mechanism by which the metastable 
configuration (trans-cisoid) decays into the stable phase (cis-transoid). In particular, we 
identify the electron-phonon configurations that mediate the nucleation process as polarons 
of the metastable phase. These polaron configurations are the equivalent of nucleation 
droplets, they correspond to saddle-points of the energy functional and the critical droplets 
(critical radius of the polarons) are identified with the transition state. The "radius" of 
the polaron (droplet) is identified as the reaction coordinate and eventually quantized as a 
collective coordinate. 
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Once these configurations are identified, we estimate the activation barriers and the decay 
rates via thermal activation and investigate the possibility that the metastable phase decays 
via quantum tunneling (quantum nucleation) at very low temperatures. Upon quantization 
of the relevant collective coordinates that describe the droplet-polaron, we provide a WKB 
estimate for the rate of quantum nucleation. Understanding the relevant electron-phonon 
configurations that mediate the metastable decay allows us to provide an estimate for the 
crossover temperature between quantum nucleation and thermal activation. 

Despite a very large body of theoretical and experimental work on quasi-one dimensional 
conductors both degenerate and non-degenerate, there has not been a consistent study of 
the fundamental aspects of the decay of the metastable phase. 

A particularly interesting and tantalizing possibility that we explore in this article is 
that the decay of the metastable phase may provide an example of mesoscopic quantum 
tunneling. 

The article is organized as follows: section II reviews the continuum model for non- 
degenerate polymers. In section III we study the constant dimerization configuration and 
establish the range of dimerization values available, scanning through values of the intrinsic 
dimerization consistent with experimental values of the gap in the stable (cis) phase. The 
energy and dimerization values obtained in this section are used to obtain the polaron 
excitations of the metastable phase. Section IV is devoted to the study of the electron- 
phonon configurations that mediate the nucleation. We construct the polaron solutions 
of the metastable phase; these are nucleating droplets, whose energy as a function of the 
"radius" of the droplet, which is identified as the "reaction coordinate", has a barrier and 
is thus identified with a saddle point in function space. The dependence of the energy 
and activation barriers upon doping (occupation of electron bound states) is analyzed in 
detail. In section V we quantize the droplet configuration via collective coordinates, obtain 
an estimate for the decay rate via thermal activation and discuss the possibility of decay 
induced by doping (electrons or holes). In section VI we study the decay of the metastable 
phase via quantum tunneling using a WKB approximation for the collective coordinates 
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of the polaron. We estimate the crossover temperature between thermal activation and 
quantum nucleation. Section VII summarizes our results and presents our conclusions and 
speculations. 



II. THE MODEL HAMILTONIAN 

The cis-transoid isomer has a slightly lower energy than the trans-cisoid isomer which 
is metastable. Because of the slightly different energy of the two configurations, there is 
an 'intrinsic' dimerization that explicitly breaks the degeneracy between the ground state 
configurations. As a result of this explicit symmetry breaking, kink excitations are not 
available. 

We consider the generalization of the SSH |J model Hamiltonian that includes the non- 
degenerate case. The discrete version of the electron-phonon Hamiltonian has been intro- 



duced and discussed by Wang, Su and Martino M, Kivelson M and Lu [TO] and is given by 



(here we assume the simpler case of a unique elastic constant for the lattice distortion) 

H = ~ £ {fa - - U 2n )](Cl, s C2n-l,s + h.C.) + 

n;s 

[t 2 - oi 2 (U 2n - U 2 n+l)](Cl n+ i )S C 2 n,s + h.c.)] \ + 



Ef(^ + l-^) 2 + ff>„ 2 (2.1) 

It proves convenient to introduce the parameters 

t = (t 1 +t 2 )/2 ; a Q = (a 1 + a 2 )/2 (2.2) 
A e = (*i - h) ; Aa = (e*i - a 2 ) (2.3) 

for the case At = 0; Aa = one recovers the original SSH Hamiltonian for trans-transoid, 
which is the degenerate isomer. Following the steps that led to the continuum version of the 
SSH model, as described by Takayama, Lin-Liu and Maki 0, one is led to the continuum 
version suggested by Brazovskii and Kirova |Q. As usual, since only electron states near 
the Fermi surface are important and the relevant phonon processes leading to the Peierls 



instability involve momenta of the order 2k f, the Fermi spectrum is linearized near the 
two Fermi points, leading to a spinor description (left and right branches). The lattice 
displacement is written as 

'A{x)\ 



Un 



"I 



4a 



(2.4) 



where the 'gap parameter' A(x) is a slowly varying function (on the scale of a lattice spacing). 
The continuum Hamiltonian, as derived by Brazovskii and Kirova ||, becomes 

K 



H 



f M 

\Z2a 2 a 



dx t ^^A 2 (x) + -^A 2 (x) 



\a 2 a 



x 



-i){v F 



aAa . . . . d . . . . . . 

A(x))cr 3 — + (A(x) + A e )a x 



(2.5) 



(2.6) 



where a is the lattice spacing a ~ 1.22A, <Ji are the Pauli matrices and ^>{x) is a spinor. 
The label s = 1, 2 corresponds to the two spin projections and plays a passive role. 

For slowly varying A(x) (on the scale of a lattice spacing) and small Aa, the term 
proportional to Aa is a small renormalization of the Fermi velocity and of the same order 
as terms that have been neglected in the derivative expansion leading to the continuum 
limit. Thus, following the arguments presented by Lu [|T0 | , we will neglect this term. The 
Hamiltonian obtained is the same as that considered by Fesser, Bishop and Campbell as 
a model for cis — (CH) X . 

Introducing the dimensionless electron-phonon coupling constant A and the bare phonon 
frequency ujq as 

K 1 



AK 



%a 2 Q a ttvfX 

the model Hamiltonian for the non-degenerate isomer becomes 



A 2 (x) A 2 (x) 



(2.7) 



H = J dx 



UqTcvfX ' hvfX 



d 



(2.8) 



We will concentrate, for the moment, on static configurations. Following Wang, Su 
and Martino || and Lu [JTOj] we will assume the same value of the parameters as for the 



trans — (CH) X degenerate proposed by Su, Schrieffer and Heeger |1J 

a = \.22A ; t = 2.5eV (2.9) 

eV eV 
« = 4.1 T ; K=21-g (2.10) 

leading to the following values of electron-phonon coupling A, band-width W, Fermi velocity 
Vp and phonon frequency 

A = 0.4077 ; W = At Q = 10 eV ; v F = 2t Q a = 6.10 eV A ; cu Q « 0.14 eV (2.11) 

in units in which h — 1. Since values of the intrinsic dimerization are not experimentally 



available (see discussions in [|8,10|) we will search for a range of values for A e such as to 
reproduce the value for the energy gap for cis — (CH) X (which is the lower energy stable 
ground state configuration); this value is approximately 2.05 eV [§,TU . 



III. CONSTANT DIMERIZATION 

The constant dimerization case corresponds to assuming A (2) = A Q (a space independent 
value). The fermions have a constant "mass" given by 

M / = A + A e (3.1) 

In the Born-Oppenheimer approximation, the electronic energy is given by the energy of the 
completely filled valence band when the fermions have the above "mass" 

E F = -2— [ kc dkJk 2 v 2 F + Mf (3.2) 

27T J -h c 

with k c , the momentum cuttoff, related to the band-width W as 

W « 2k c v F 

Incorporating the elastic energy, we obtain the energy per site 



E 



A 



2 r w/2 , : 

°- - 2 / dzJz 2 + M 2 
Jo 3 



The extrema of this energy function are determined by the "gap equation'' 



(3.3) 



A 2 
A 



W/2 

= / dz 



M 



f 



z 2 + M 2 



(3.4) 



For the values of the parameters given by Su, Schrieffer and Heeger for the degenerate case 
(trans), we show in Figure (1) the energy per site (in eV) as a function of A (in eV) for the 
representative value A e = 0.02 eV . As the value of A e is increased, the stable minimum 
becomes deeper and the metastable minimum becomes shallower and closer to the origin. 
When A e > 0.07(eV) the metastable minimum disappears altogether. Numerically, the 
optimum range of A e that predicts a gap for the continuum theory with approximately the 
same error as the continuum theory prediction for the degenerate case, with the SSH values 
for the parameters (about 10 %) [|1^], is 0.02 < A e < 0.07. In this range we find: 



A e (eV0 


A _(eV) 




A 0+ (eV) 




0.02 


-0.835 


-1.610 


0.896 


-1.621 


0.04 


-0.800 


-1.604 


0.923 


-1.626 


0.06 


-0.760 


-1.596 


0.947 


-1.632 



(3.5) 



Where A c ± correspond to the values of A Q at the metastable (-) and stable (+) minima 
with energy per site E±/N respectively. 

Thus, the value of the gap predicted from the global minimum (2A 0+ ) for the stable 
configuration is fairly close to the observed gap 2.05 eV for the lower energy cis — (CH) X 
configuration (again the discrepancy is of the same order as the discrepancy between the gap 
predicted by the continuum model and the observed value in the degenerate case (trans), with 



the value of the parameters chosen by SSH [0). For values of 0.02(eV) < A e < 0.07(eV^) 
curves very similar to figure (1) are obtained for the energy per site. 
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IV. POLARON SOLUTIONS 



The non-degenerate case does not permit soliton (kink) type excitations (except if a 
very special theoretical possibility is fulfilled but allows the possibility of polaron 

(bag) excitations. Polarons are topologically trivial electron-phonon field configurations 
that correspond to distortions of the lattice in which there are trapped electrons (bound 
states). Because of charge conjugation symmetry, bound states must appear symmetrically 
with respect to the middle of the gap. 

For polaron solutions, the phonon field profile reaches asymptotically the values of the 
minima of the effective energy for the constant dimerization case, that is the minima in figure 
1. The polaron solutions corresponding to the global (stable) minimum, had already been 
studied in references( |5i[7|, [10| , [l5| ), these solutions reach A Q+ asymptotically when \x\ — > oo. 
These solutions are parametrized by two collective coordinates: the center of mass of the 
polaron (position), associated with translational invariance, and the radius. For the polaron 
solutions in the stable (lower energy) phase, the polaron energy is a linearly increasing 
function of the radius, for large radius. These polarons are "confined". The reason is 
simple: for large radius, the polaron is exploring a region in function space that is very 
close to the higher (metastable) minimum, thus increasing the "volume" energy as a linear 
function of the radius. 

In this article, we are interested in the polaron configurations in the metastable phase. 
We argue that these are the relevant electron-phonon configurations responsible for the decay 
of the metastable state. 

The argument is as follows: let us assume that a polaron solution exists in the metastable 
phase. This configuration will reach the values of the metastable minimum for the constant 
dimerization case (A _) asymptotically as \x\ — > oo. Its profile again will be that of a "bag", 
inside which the field configuration will sample the lower energy minimum in a region in 
space given by the radius of the polaron. However, the stable minimum has lower energy and 
the system gains volume energy by increasing the radius. There is a price in elastic energy 
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determined by the gradients of the field configuration, but in one dimension this "surface 
energy" is independent of the radius for large radius and bound by a value close to twice the 
mass of the kink-antikink pair that is the asymptotic state of the polaron for infinite radius. 
For small radius, the elastic term wins out and the energy is an increasing function of the 
radius; for large radius, the elastic energy will saturate and the volume energy will dominate; 
in this case the energy is a decreasing function of the radius. Thus, if a polaron solution 
exists with asymptotic values for the dimerization in the metastable phase, there must be a 
critical radius. This configuration is equivalent to the critical "droplet" configuration that 
Langer introduced to explain the statistical decay of a metastable state in first order phase 



transitions [17j and the "radius" of the droplet is identified with the reaction coordinate. 

We will argue below that these polaron configurations will be saddle points of the energy 
functionals as they will be characterized by an unstable mode, a zero mode and fluctuations 
with positive frequencies. 

To actually construct the polaron configurations, we could invoke the large body of 



work on the exact integrability of the continuum theory [|5|, [TD|JTT| , [T5H r5|] , but we prefer to 
go through some of the details of the calculation, as there are some subtle but important 
features that must be addressed. 

Borrowing from the known results for the degenerate case fll| and for the stable minimum 
|5|j7|JT0[| (cis) we propose the static polaron solution as 

A p (x; x cm ; x Q ) = A - K Q v f {tanh [K Q (x - x cm + x Q )} - tanh [K a (x - x cm - x Q )]} (4.1) 



The fermionic potential becomes reflectionless for the integrability condition JTT| , P^ r5 

tanh[2if x ] =^ (4.2) 
Mf 

M f = A + A e (4.3) 

In the expressions above, A D can be either one of the minima of the energy functional 
for the constant dimerization case, i.e., A G ±. 

Thus, the polaron profile is parametrized in terms of x Q and x cm , the "center of mass"; 
later we will treat these as "collective coordinates". The "radius" of the polaron is 2|x |. The 
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center of mass coordinate x cm reflects the underlying translational invariance, the energy is 
independent of this coordinate. 

Before proceeding further to the computation of the polaron energy, it proves illuminating 
to understand some features of the polaron solution and the integrability condition: i) the 
polaron solution (and consequently the electronic spectrum) is an even function of K . 
We choose K a > 0; ii) for a given x Q , the integrability condition (|4.2|) determines K a ; this 
condition also determines that the sign of x Q is given by the sign of Mf. This is an important 
point. The stable minimum (A n +) always corresponds to Mf > 0, and consequently x a > 0. 
The polaron of the stable phase decreases the value of the local dimerization in the region 
of the lattice distortion, thus sampling a region of higher energy density (trans). When x Q 
becomes large the energy will grow linearly with x . This is the "confining" mechanism 
PJ7|JT0[| in the stable phase. 

Now consider the case in which A G is the value at the metastable minimum (A c _). In this 
case Mf < and consequently x < 0. The polaron in this phase increases the value of the 
local dimerization in the distortion region, thus sampling regions with lower energy density. 
Thus for large | x Q | the energy will decrease linearly with | x | and it becomes favorable for 
the system to produce large polarons. But, clearly, there will be an energy barrier to do so 
because, for small radius, the elastic energy will gain as will be shown below. 

iii) The integrability condition (|4.2| ) yields non-trivial solutions for K only when 2\x D \ > 
v f /\Mf\. 

For 2\x a \ < Vf/\Mf \ there is no polaron solution and the only available solution is that 
for constant dimerization. The physical reason for this is that Vf/\Mf \ is the Compton 
wavelength of the fermions, and there cannot be bound states localized within a region in 
space smaller than the Compton wavelength. 

In the other limit when \x Q \ — > oo, K Q — > \Mf\/vf, the electronic bound states (see below) 
merge at the middle of the gap (zero energy). In this case, the polaron looks like a widely 
separated kink-antikink pair, each with a localized fermionic zero mode; the two bound states 
correspond to the symmetric and antisymmetric combinations of the fermionic "zero mode" 
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wave functions, split-off in energy because of the (small) overlap of their wave-functions. 

It remains to compute the energy of the polaron. There are two contributions: the elastic 
contribution from the term A 2 (x) in the Hamiltonian, and the electronic energy. Due to 
the reflect ionless nature of the potential, the electronic spectrum is known exactly P,|TT 



there are two electron bound states at energies ±co = ±yM| — K%Vf (the wave functions 
of these bound states are localized at x cm ± x ) and the conduction and valence continuum 



states with dispersion E c = ±yk 2 vj + Mj. The valence and conduction bands are both 
depleted by one state per spin degree of freedom, the scattering phase shifts of the valence 



band states are 11 



S(k) = 2 tan" 1 



Ko 
k 



(4.4) 



The form of the bound states and continuum electronic wave functions has been given exactly 
by Campbell and Bishop JT1J for the polarons in the degenerate case. Their results can be 



applied vis a vis to our case because they are properties of the electronic "Dirac" equation 
in the presence of a spatially varying background configuration whose profile is the same as 
in the degenerate case. The reader is referred to those articles for more details. 

The difference in energies between the polaron and the constant dimerization case is 
written as 

6E = 8E l + 6E el (4.5) 
SEt = — - J dx{A 2 p (x; x cm - x ) - A 2 } (4.6) 
8E d = (n+ - n-)u - 2 $> P (A;) - co°(k)) (4.7) 

k 

where we have assumed that the positive and negative energy bound states are occupied 
with n + , n~ electrons, respectively (n^ = 0,1,2), the factor 2 accounts for the two spin 
projections and u p, ° are the energies of the states in the valence band, with and without the 
polaron, respectively. We find 

= 8A,A> -4/^, 
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4 4 f A c 

5E e i = (n + -n )u + -K Q v f + -K D v f — - 

IX 7X \ AM 

K o v f 

UJ 



S. 



4 (TX 

+ —u - - tan 

7T V Z 



(4.9) 



The first term in ( |4.8f ) has the important physics that we we were looking for. Whereas 
the rest of the terms in fl4.8|J4.9|) reach a constant value when \x \ — > oo, it is this first term 
that dominates the contribution to the energy in this limit. 

As argued above, for the polaron in the metastable state x Q < 0, whereas for the polaron 
in the stable state x Q > 0. Thus, in the metastable state, as the size of the polaron becomes 
large, the energy becomes large and negative, linearly with the size. In the stable state it 
becomes large and positive linearly (this is the confining mechanism found in [|5|,fr|,|T0f ) . 

To compare the energy in the metastable state with the results quoted in the literature 
for the stable case, it proves convenient to introduce the variables 9 and 7 as 

K v f = \M f \ sm{6) ; u =\M f \ cos(fl) ; < 9 < ^ (4.10) 

Finally, the energy difference in terms of these variables is given by 

5E=-\M f \ |7tanh" 1 [sin(0)] + sin(0) - 7sin(0) + ^(n + - n~) cos(#) + 

oos(6l) (I -6l)}. (4.12) 

It is important to recognize that, for the minimum at the metastable state, 7 < 0. 

This expression is similar in form to the one found in references ( ]5|,[7 , |10 , [15 1 ) with the 



only difference being the subtleties associated with the signs for the metastable case. 

To contrast the polaron solutions in the metastable state to those in the stable state, let us 
look at the extremum condition obtained from the energy function( |4.12| ). The integrability 
condition ( |4.2| ) relates K Q (and u ) and phase shifts to x a . The value of x a is obtained from 
the extremization condition of the energy. Because x Q is a monotonically increasing function 
of 9 in the interval < 9 < tt/2, it proves more convenient to extremize with respect to 9. 
We obtain the equation 
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0-|7|tau(0) = -(n + -n_+2). (4.13) 

For I7I > 1 the only solution is the trivial one, 8 = 0, however this value of I7I is not within 
the allowed range of parameters that describe a metastable situation. For I7I < 1 there are 
several interesting possibilities: 

i) for (n+ — nJ) = —2 there are always two solutions, one corresponding to the trivial 
case 9 = and another non-trivial solution. This is in marked contrast with the stable 
situation for which, at these values of the occupation for the bound states, there is only 
the trivial solution. Thus, in the metastable case there is always a polaron solution for the 
ground state configuration n + = 0,n_ = —2. Because the energy difference between the 
constant dimerization and the spatially varying case vanishes at 9 = and grows linearly 
with negative slope for large \x Q \ (6 — > 7r/2), this polaron solution corresponds to a maximum 
of ( |4.12|) as will be shown explicitly below. 

ii) When n + — n_ + 2 > non-trivial solutions are available for a particular range of 
parameters and will be studied numerically below. However, we find that whenever non- 
trivial solutions are available, they always appear in pairs. A similar analysis as presented 
above reveals that the solution with the smallest value of 9 corresponds to a minimum and the 
largest to a maximum of ( |4.12j ). These solutions correspond to polaron-like configurations 
in which a lattice distortion traps electrons in positive and negative energy bound states 
and are the analog of the polaron solutions found in references ( [5,[7|,|IDJ) for the stable 
configuration. 

Figures (2.a,b,c) show the energy as a function of r Q = MfX Q /vf for n + — n_ = —2; — 1; 
for A = 0.4077, A e = 0.02 eV . For n + — n_ = 0, there is no minimum or maximum, i.e. 
there are no solutions to the extremum equation ( |4.13|) ; we will refer to this case as the 
"dissociation curve", for reasons that will become clear later. Similar curves are obtained 
for A e = 0.04 eV , where again there are polaronic solutions only for n + — n_ = —2; —1 
and dissociation curves for n + — ri_ > . For A e = 0.06 eV only for the ground state 
configuration n + — n_ = —2 there is an extremum "polaron" solution corresponding to a 



14 



maximum of the energy, for all other values of n + — n_ we find dissociation curves without 
extrema. From the figures, one finds that in all cases the typical size of the polaron solutions 
is2|a; | « 20 - 40A. 



V. DECAY OF THE METASTABLE STATE: THERMAL ACTIVATION 

We are now in condition to study the mechanism that leads to the decay of the metastable 
state. The polaron configurations corresponding to the maximum of the energy (|4.12 ) are 



identified as Langer's critical droplets OS] or the "transition state" p9| . These con- 
figurations correspond to a saddle point of the energy functional in the multidimensional 
configuration space. The "critical droplet" is determined by the value x a = x* at which the 



energy Q4.12|) is a maximum (the maxima in figures 2.a-c). If we consider small fluctuations 



around the critical droplet (transition state) configuration, we find that there is one "zero 
mode" corresponding to translations, one unstable mode corresponding to dilation of the 
radius of the droplet and, presumably, infinitely many perpendicular directions with positive 
real frequencies, for small oscillations around the droplet configuration. The presence of the 
translational and unstable mode is easy to understand. Because of translational invariance 
the energy does not depend on the position of the polaron, that is, the coordinate x cm is 
cyclic. Then, a fluctuation around the polaron solution of the form 

<W = a a Qx (5.1) 
with a Q constant and small, corresponds to a shift of the polaron position 

but, by translational invariance, the energy functional is invariant under such a shift. Thus, 
the function 

dA p (x) , . 

J [X) X cmj X J yo.o) 
OX 
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and the collective coordinate a Q determine a "flat" direction in functional space associated 
with the "zero mode". Notice that the function f Q (x) is antisymmetric and has one node. 

Now consider a small fluctuation of the polaron solution around the value of x a = x* 
with 2x* the value of the radius of the polaron corresponding to the maximum of the energy 
function (4.12). This fluctuation is determined by the function 



S ^ x ) = Q i £ Uo=x% (5.4) 



dx 



o 



For small a\ 

A p (x; x cm ; x* Q ) + 8i(x) w A p (x; x cm ; x* a + a x ) (5.5) 

Since x* Q corresponds to a maximum of the energy functional, the coordinate a\ determines 
(locally) an unstable direction in function space around the extremum solution. It is im- 
portant to notice that the function j\ (x) is symmetric (nodeless) and is thus orthogonal to 

fo(x). 

Generalizing Holstein's |20|-p2]j treatment of the large polaron to incorporate Xq as a 



collective coordinate, we write the quantum expansion around the polaron solution as 

A(x) = A p (x - x cm ; x ) + aifi(x ~ x cm ; x Q ) (5.6) 

i>i 

where the functions fi{x) are constrained to be orthogonal to f (x) and fi(x) and correspond 
to the stable modes of oscillations around the polaron solutions. 
From the Hamiltonian (12.81) we find 



with 



M 2 M 2 

H = -f±\ m + -^x\ + E( Xo ) + ■■■ (5.7) 



= JdxjdAp/dx] 2 
\nv f u; 2 Q 

= Jdx[dAp/dx f 

\TTVfU 2 Q 
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and where the dots represent the coordinates associated with other functional directions 
and interactions, and E(x ) is the Born-Oppenheimer energy (|4.12|) in terms of 9(x ). The 
masses M Q and Mi are numerically found to have a very weak dependence on x D for values 
of x in the region of interest (near the maxima and minima). We find 

M « Mi « 4.2( eV A 2 )" 1 (5.10) 



Following Langer |17j and Affleck [18|] (see also |19| ) the decay rate of the metastable 
state is obtained as 

„ n 



2vr" 



-ImF 



(5.11) 



where Im F is the imaginary part of the analytically continued free energy computed in the 
saddle point approximation around the polaron solution corresponding to the maximum of 
the energy function (i.e. the transition state) and Q is the (unstable) frequency at the top of 
the barrier along the functional direction determined by dilation of the radius of the polaron. 
Near the maximum of the energy function, along the unstable direction, the Hamiltonian 
becomes 



H 



E ( x o) + v 1 ™ + ~ir x o - -tt^ ( X o - x ) 



(5.12) 



2 cm 2 o 2 

where E(x*) is the energy at the top of the barrier (the maximum) and the dots again stand 
for the stable modes and possible interactions. We find numerically the following values 
for the energies at the maxima (E(x*)), minima (E min ), bound state energies u ± at the 
maxima (+) and minima (— ) respectively, and unstable frequency Q at the top of the barrier 
(An = n + — n_). 



A e (eV) 


An 


E(x* )(eV) 


E min (eV) 


u 0+ (eV) 


Lu -(eV) 


Q(eV) 


0.02 


-2 


0.792 




0.032 




0.031 


0.02 


-1 


0.838 


0.722 


0.070 


0.535 


0.103 


0.04 


-2 


0.573 




0.066 




0.077 


0.04 


-1 


0.672 


0.655 


0.171 


0.425 


0.070 


0.06 


-2 


0.393 




0.102 




0.149 



(5.13) 
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Thus, we see that the unstable frequencies at the top of the barriers are typically of the 
same order of magnitude as the bare phonon frequencies. 

The decay rate is thus obtained by calculating the partition function at the saddle point, 
the unstable mode is treated via an analytic continuation, normalized to the partition func- 
tion in the metastable well. The zero mode leads to a volume dependence (L) [JT7]] and we 
find 



1 [M k B Tf 2 e W 



[K]. (5.14) 



L 2?r sin(^ /k B T) 

[K] is the ratio of the square root of the determinants for the stable modes around the 
extremum polaron solution and the constant dimerization solution for the metastable well at 
finite temperature. The computation of this ratio of determinants is beyond our capabilities, 
it is a dimensionless number, presumably of order one, as the frequency scales in both cases 
are of the same order of magnitude. The singularities in ( |5.14|) at T = Q/2nk B n have been 



discussed by Wolynes J^J and the reader is referred to that article for details. In our case for 
n^O these values of the temperature are below the crossover value for thermal activation 
to quantum nucleation (see discussion below) and the rate given by ( |5.14| ) is no longer 



applicable. As shown in the table above, OAeV < E(x* Q ) < 0.8eV thus for temperatures 
T « 100 — 200K the Kramers- Arrenhius factor in ( |5.14| ) is fairly large and the lowest energy 
state in the metastable phase (trans), the constant dimerization, is fairly long-lived. 



A. Induced decay 

The lowest energy configuration in the metastable state is that for constant dimeriza- 
tion with the valence band completely filled. Upon doping with electrons (or holes) the 
metastable phase, the energetically most favorable configuration corresponds to the forma- 
tion of a polaron in the Born-Oppenheimer surface corresponding to n + — n + = — 1 with 
either one electron in the bound state at +u a and two electrons in the bound state at — uj 
(electron polaron) or one electron in the bound state at — uo (hole polaron). On this energy 



18 



surface, the lowest energy polaron configuration corresponds to the minimum (available for 
A e < 0.06 eV ). 

From the table above, we see that the difference in energy between the minimum and 
maximum polaron configuration on these surfaces (E(x*) — E m i n ) is typically an order of 
magnitude smaller. In this situation, the activation energy barriers are of order AE « 
0.02 — 0.1 eV , and the decay rate is dramatically enhanced. 

Upon further doping, the Born-Oppenheimer surface may change to n + — > and 
the dissociation curve is reached for almost all phenomenologically available values of A e . 
In particular, this is always the case for bipolarons n + = n_ for the range of parameters 
consistent with the phenomenology of cis-polyacetylene; this is another major difference 
with the polarons in the stable phase. In this case, there are no barriers and the system 
decays instantaneously. The decay process, in this situation, cannot be studied under the 
assumptions of quasi-equilibrium implied in the treatment of metastable decay and a time 
dependent non-equilibrium treatment would have to be used. 

The first excited and dissociation curves may also be reached by absorption of photons 
of energy %v = \Mf\ + u and electronic transitions to the bound state with energy +u . 
An absorption peak in this energy range will be the telltale of electron bound states, just as 
absorption peaks in the degenerate case reveal the existence of electrons bound to solitons 
f24| . Since uj is typically very small, these states may appear as mid-gap states at %v ~ leV. 
However, there is no room for confusion with bound states on solitons, because solitons are 
not available in the non-degenerate case. 

VI. DECAY VIA QUANTUM TUNNELING 

At very low temperatures (to be quantified later), we expect that the metastable state 
will decay via quantum tunneling. In order to understand quantum tunneling in multidi- 
mensional space, one must search for the configurations that extremize the classical action 
in imaginary (euclidean) time. These configurations constitute the "most probable escape 
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path" [21,0 , and are solutions to the classical equations of motion in euclidean time [ 271 , 



dubbed "bounces" [28-32 



In terms of the collective coordinate x Q representing the radius of the polaron, this "most 
probable escape path" or "bounce" corresponds to the classical trajectory in the inverted 
potential —E(x ) between the two classical turning points corresponding to the energy of 
the metastable state. For the ground state Born-Oppenheimer curves (n + — ri- = —2) this 
energy is zero, whereas for the higher energy surfaces (n + — = —1) the energy is that of 
the minimum and marked with a dashed line in figure (2.b). This approach was previously 



used by Kivelson and Sethna [Q to study photoinduced soliton-pair production. 

To exponential accuracy in the semiclassical WKB approximation, the decay rate per 
unit length is given by 

r 



— m exp 
L 



■4 



(6.1) 



with 5*0 the action of the classical trajectory in euclidean time at energy E between the two 
classical turning points x_; x + : 



So = / dx J2M [E(x ) - E] (6.2) 



Again, the total decay rate will have prefactors in front of the WKB exponential in fl6.1|) ; 
this prefactor involves the ratio of the determinants of the quadratic fluctuations. We are 
unable at the moment to calculate the prefactor and content ourselves with an estimate of 
the exponential, which is the leading term in the semiclassical approximation. 

We can provide a rough estimate of the crossover temperature between thermal activation 
and quantum nucleation by comparing the WKB factor to the Arrhenius- Kramers activation 
factor, thus obtaining the approximate estimate for the crossover temperature 

hAE 



2S n 



(6.3) 



where AE is the activation barrier (E(x ) for the ground state or E(x Q ) — E min for the 
higher energy surface). The action S Q is calculated numerically and we find: 
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A e (eTO 


An 


So 

h 


k B T c (eV) 


0.02 


-2 


102.7 


0.004 


0.02 


-1 


11.3 


0.005 


0.04 


-2 


41.3 


0.007 


0.04 


-1 


1.4 


0.006 


0.06 


-2 


21.0 


0.009 



(6.4) 



Thus, we see that, typically, the crossover temperature is of the order of T c pa 40 — 100K 
depending on the value of the intrinsic dimerization, and may be amenable to experimental 
realization. 

Another (but related) criterion for the crossover temperature is obtained from the 
quasiequilibrium approach p2 |. It is based on the properties of the "bounce" solution in 
Euclidean time corresponding to a classical trajectory with period (3h in the inverted poten- 
tial well. These trajectories give the leading semiclassical contribution to the equilibrium 
partition function. The smallest period for trajectories in the inverted potential corresponds 
to the harmonic oscillations with frequency Q at the bottom of the (inverted) well (see 



equation |5.12| ). Thus the maximum temperature for which there are "bounce" solutions in 
Euclidean time is given by 



(6.5) 



leading to an estimate for the crossover temperature jTJ 

TiVL 



2nk 



B 



(6.6) 



With the values of Vt from ( |5.13| ), we find the estimate for crossover temperature based 
on this criterion to be 50f^ < T c < 100K which is consistent with the previous estimate 
based on the comparison of the Arrenhius-Kramers and WKB factors. For T < T c the 
decay process will be dominated by quantum nucleation and will be fairly insensitive to 
temperature, whereas for T > T c the decay will be dominated by thermal activation and 
the rate will depend strongly on temperature via the Arrenhius-Kramers factor in the decay 



21 



rate. Again, quantum decay will be accompanied by the presence of electron bound states 
at energies ±u and may be detected via absorption peaks, just as in the degenerate case, 
where the absorption peak at the middle of the gap signals the presence of solitons. 

Quantum nucleation can also be induced by doping, thus moving from one Born- 
Oppenheimer surface, for example the ground state, to the next for which n + — n_ = — 1. 
We see that the height and width of the barrier for this surface are much smaller, resulting 
in a ten-fold reduction of the WKB action. Thus, the transmission probability, i.e. the 
nucleation rate, is greatly enhanced. Upon further doping, the dissociation curve is reached 
and there is no quasi-stationary state; one must resort to a full real time non-equilibrium 
calculation of the decay probability in this case. 

The decay through quantum nucleation occurs via the spontaneous production of a fairly 
large droplet; this is a tunneling process. The size of the droplet will be roughly determined 
by the value of the coordinate \x a \ at the classically forbidden turning point below the 
barrier, i.e. the "escape point" . For the lowest energy surface, this value is typically w 40A 
whereas for the first excited surface it is about ~ 20A. Because of these large values of the 
radius of the quantum droplet, the bound states will be very close to the middle of the gap 
and this quantum droplet looks like a kink-antikink pair widely separated. 

The telltale signal for the decay of the metastable (trans-cisoid) state, either via thermal 
activation or quantum nucleation, will be electron bound states with energies very near 
the middle of the gap (because for the typical sizes of the polarons, the energy oj is very 
small) and detected through absorption peaks at about %v « leV. These peaks could not be 
confused with solitons as these topological excitations are not available in the non-degenerate 
isomers. Furthermore, the different regimes for thermal activation and quantum nucleation 
may be separated by plotting the logarithm of the nucleation rate versus 1/T: activated 
processes lead to an approximate straight line with negative slope, whereas for quantum 
nucleation there should be a flat plateau. Our analysis for the crossover temperature leads 
to the prediction that such a plot should have a plateau for temperatures smaller than about 
AOK. 
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We must say that at the moment we are not aware of experiments that have either looked 
at or reported on these possibilities. 

VII. CONCLUSIONS AND SPECULATIONS 

The motivation of this work was to study the mechanism responsible for the decay of the 
metastable (trans-cisoid) configuration to the stable (cis-transoid) isomer in non-degenerate 
polymers. We have identified the electron-phonon configurations that play the role of the 
nucleation droplets as polarons in the metastable phase, and have studied their properties 
within the range of parameters for the intrinsic dimerization compatible with a gap of about 
2 eV in the electronic spectrum for the cis-isomer. In the undoped, lowest energy metastable 
configuration, typical barrier energies are of the order of about 0.4 — 0.8 eV depending on the 
intrinsic dimerization parameters. Upon doping with electrons or holes, higher energy Born- 
Oppenheimer surfaces become available and activation barriers become smaller by almost 
an order of magnitude; thus, upon doping, the metastable decay is induced and the rate 
is enhanced dramatically. For large doping, dissociation curves become available without 
activation barriers; in this case a full time dependent non-equilibrium study will be needed 
to understand the decay. 

We obtained an estimate for the decay rate via thermal activation by treating the relevant 
coordinates of the polaron (center of mass and radius) as collective coordinates and evaluated 
numerically the necessary activation barriers and unstable frequencies to obtain the rate 
approximately. We also performed a semiclassical WKB calculation to estimate the decay 
rate via quantum nucleation and to establish the crossover temperature that separates the 
regimes between thermal activation and quantum nucleation. We found these temperatures 
to be approximately T c < 40 K. 

We speculate that metastable isomers in non-degenerate polymers may be candidates to 
study mesoscopic quantum tunneling. Typical values for the sizes of critical droplets (po- 
larons) are 20 — 40A The telltale signal that the metastable decay is mediated by polarons 
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will be electron bound states with energies very near the middle of the gap. These states 
may be detected via absorption peaks at about 1 eV, just as the electron states bound to 
solitons in the degenerate case. There is no possibility to confuse the metastable polarons 
with solitons as the latter are not available excitations in the non-degenerate case. Further- 
more, our rough estimate for the crossover temperature should be amenable to experimental 
confirmation by plotting the logarithm of the nucleation rate versus 1/T. In the thermally 
activated regime this should be almost a straight line with negative slope, the curve should 
flatten and reach a plateau at temperatures T < 40K, regime in which the nucleation rate 
is insensitive to temperature signaling quantum nucleation. 
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Figure Captions: 

Figure 1: (eV) vs. A (eV) for A = 0.4077; A e = 0.02 eV . 

Figure 2. a: SE(r )(eV) (equation (|4.12|) ) vs. r = MfX a /vf, for n + — n_ = —2. The 
"bounce" trajectory has zero energy. 



Figure 2.b: 5E(r )(eV) (equation ( |4.12|) ) vs. r Q = MfX Q /vf, for n + — n_ = — 1. The 



dashed line represents the energy for the "bounce" trajectory. 

Figure 2.c: 5E(r )(eV) (equation ( [4.12|) ) vs. r a = MfX /vf, for n + — n_ = —2; — 1; 0. 
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